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LONG-RANGE  PROPAGATION 


LONG-TERM  GOAL 

Our  long-term  goal  is  to  provide  a  more  complete  understanding  of  the  forward  propagation  of 
acoustic  pulses  in  range  dependent  deep  ocean  environments  at  multi-megameter  range  scales. 

SCIENTIFIC  OBJECTIVES 

There  are  two  primary  scientific  objectives  of  this  work:  1)  to  quantify  the  relevant  ocean 
environment  parameters  that  determine  the  range  where  wavefield  predictability  is  lost  (for  a 
given  transmitter  center  frequency  and  bandwidth),  and  2)  to  develop  a  model  based  on 
geometric  acoustics  that  facilitates  our  understanding  of  long  range  (multi-megameter)  acoustic 
propagation.  Both  of  these  objectives  play  a  role  in  the  ultimate  goal  of  understanding  the  limits 
of  'ray-based'  acoustic  tomography. 

APPROACH 

Traditionally,  long  range  propagation  of  sound  in  the  range  dependent  ocean  waveguide  has  been 
treated  using  concepts  developed  from  the  field  of  waves  in  random  media.  Our  approach 
attempts  to  use  these  traditional  treatments  as  well  as  geometric  concepts  from  the  dynamical 
systems  field  to  give  insight  as  to  what  to  expect  upon  propagating  a  single  acoustic  pulse,  where 
the  ocean  environment  is  taken  to  be  deterministic.  For  sound  speed  fluctuations  due  to 
mesoscale  structure,  it  is  anticipated  that  a  deterministic-based  treatment  will  be  realizable  in  the 
near  future. 

The  ray  equations  form  a  dynamical  system  that  is  known  to  be  nonintegrable  —  individual  ray 
trajectories  exhibit  exponential  sensitivity  to  the  environment  —  when  the  ocean  has  range- 
dependence.  At  any  given  range,  a  family  of  rays  emanating  from  a  point  source  form  a 
Lagrangian  manifold  in  a  phase  space  that  is  composed  of  each  ray  trajectory's  position,  angle, 
and  travel  time.  Folds  of  the  manifold  in  this  phase  space  dictate  where  caustics  reside.  Our 
numerical  approach  consists  of  solving  the  three  dimensional  ray  dynamical  system  through  a 
quasi-realistic  ocean  environment  composed  from  the  output  of  a  eddy  resolving  ocean  global 
circulation  model.  The  locations  where  folds  in  the  Lagrangian  manifold  occur  are  to  be 
recorded,  thus  indicating  the  range  dependence  of  both  horizontal  and  vertical  caustics  generated 
due  to  the  nonintegrable  nature  of  the  dynamical  system.  This  approach  will  provide  insight  as 
to  the  importance  of  horizontal  refraction  due  to  mesoscale.  The  model  is  designed  to  be 
extensible,  so  that  smaller  scale  ocean  phenomena  (e.g.  internal  waves)  can  be  included. 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  1997  2- REPORT  TYPE 

3.  DATES  COVERED 

00-00-1997  to  00-00-1997 

4.  TITLE  AND  SUBTITLE 

Investigation  of  Acoustic  Wavefield  Dynamics  in  Quasi-realistic  Ocean 
Environments 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Pennsylvania  State  University, Earth  System  Science  Center, University 
Park, PA, 16802 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

3 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


WORK  COMPLETED 

We  have  completed  an  investigation  based  on  deriving  a  two  dimensional  horizontal  plane  wave 
parabolic  equation  that  follows  from  the  adiabatic  mode  approximation.  The  corresponding  ray 
dynamical  system  is  analyzed  both  theoretically  and  numerically  for  idealized  mesoscale 
environments.  This  work  is  guiding  us  in  the  three  dimensional  numerical  modelling  effort  as  to 
the  importance  of  horizontal  refraction,  at  least  for  the  near-axial  rays.  A  paper  detailing  this 
effort  has  been  submitted  to  JASA. 

Using  a  semi-classical  extension  of  this  same  model,  and  a  horizontal  parabolic  equation  model, 
extensive  numerical  simulations  were  performed  to  investigate  the  breakdown  of  ray  theory. 
Direct  wavefield  comparisons  indicate  excellent  agreement  out  to  ranges  of  four  megameters. 
Limits  of  predicting  the  wavefield  beyond  this  range,  are  seen  to  be  constrained  by 
computational  difficulties  regarding  the  resolution  of  the  Lagrangian  manifold,  as  opposed  to  a 
global  breakdown  of  ray  theory  due  to  the  accumulation  of  phase  errors.  The  preliminary  results 
of  this  work  have  been  presented  at  the  June  1997  ASA  meeting  in  State  College,  PA. 

In  regards  to  the  vertical  plane,  a  numerical  method  has  been  developed,  based  on  the  adiabatic 
approximation,  for  determining  the  range  dependence  of  vertical  scattering  due  to  internal  wave 
induced  sound  speed  fluctuations.  A  simplified  version  of  this  method  has  been  described  in  a 
paper  submitted  to  JASA. 

A  three  dimensional  ray  model  is  currently  under  development  on  the  Cray  T3E  supercomputer 
at  the  Arctic  Region  Supercomputing  Center  at  the  University  of  Alaska.  At  this  time,  the  model 
is  still  in  the  testing  phase,  and  our  first  results  quantifying  the  importance  of  horizontal  vs 
vertical  multiple  forward  scattering  will  be  presented  at  the  June  1998  ASA  meeting  in  Seattle, 
WA. 

IMPACT/APPLICATION 

The  three  dimensional  ray  code  is  expected  to  contribute  to  the  understanding  of  a  planned  long 
range,  deep  water,  billboard  array  experiment,  to  be  conducted  by  P.  Worcester  and 
collaborators.  Our  results  concerning  wavefield  predictability  are  anticipated  to  have  a  strong 
impact  on  ray-based  tomographic  methods. 

RELATED  PROJECTS 

Our  work  involves  extensive  collaborations  with  the  following  individuals:  M.  Brown,  F. 
Tappert,  and  S.  Tomsovic. 
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